X 



>ni 6 



BAR4R-PUB-08/017 
SLAC-PUB-13305 

Measurement of Ratios of Branching Fractions and CP- Violating Asymmetries of 

B± -> D*K± Decays 

B. Aubert, 1 M. Bona, 1 Y. Karyotakis, 1 J. P. Lees, 1 V. Poireau, 1 E. Prencipe, 1 X. Prudent, 1 V. Tisserand, 1 
J. Garra Tico, 2 E. Grauges, 2 L. Lopez afc , 3 A. Palano a6 , 3 M. Pappagallo a6 , 3 G. Eigen, 4 B. Stugu, 4 L. Sun, 4 
G. S. Abrams, 5 M. Battaglia, 5 D. N. Brown, 5 R. N. Cahn, 5 R. G. Jacobsen, 5 L. T. Kerth, 5 Yu. G. Kolomensky 
G. Kukartsev, 5 G. Lynch, 5 I. L. Osipenkov, 5 M. T. Ronan, 5 'B K. Tackmann, 5 T. Tanabe, 5 C. M. Hawkes, 6 N. Soi 
A. T. Watson, 6 H. Koch, 7 T. Schroeder, 7 D. Walker, 8 D. J. Asgeirsson, 9 T. Cuhadar-Donszelmann, 9 B. G. Fulsoi 
C. Hearty, 9 T. S. Mattison, 9 J. A. McKenna, 9 M. Barrett, 10 A. Khan, 10 L. Teodorescu, 10 V. E. Blinov, 11 
A. D. Bukin, 11 A. R. Buzykaev, 11 V. P. Druzhinin, 11 V. B. Golubev, 11 A. P. Onuchin, 11 S. I. Serednyakov, 11 
Yu. I. Skovpen, 11 E. P. Solodov, 11 K. Yu. Todyshev, 11 M. Bondioli, 12 S. Curry, 12 I. Eschrich, 12 D. Kirkby, 12 
A. J. Lankford, 12 P. Lund, 12 M. Mandelkern, 12 E. C. Martin, 12 D. P. Stoker, 12 S. Abachi, 13 C. Buchanan, 13 
J. W. Gary, 14 F. Liu, 14 O. Long, 14 B. C. Shen, 14 <B G. M. Vitug, 14 Z. Yasin, 14 L. Zhang, 14 V. Sharma, 15 
C. Campagnari, 16 T. M. Hong, 16 D. Kovalskyi, 16 M. A. Mazur, 16 J. D. Richman, 16 T. W. Beck, 17 A. M. Eisner, 17 

C. J. Flacco, 17 C. A. Heusch, 17 J. Kroseberg, 17 W. S. Lockman, 17 T. Schalk, 17 B. A. Schumm, 17 A. Seiden, 17 
L. Wang, 17 M. G. Wilson, 17 L. O. Winstrom, 17 C. H. Cheng, 18 D. A. Doll, 18 B. Echenard, 18 F. Fang, 18 

D. G. Hitlin, 18 I. Narsky, 18 T. Piatenko, 18 F. C. Porter, 18 R. Andreassen, 19 G. Mancinelli, 19 B. T. Meadows, 19 
K. Mishra, 19 M. D. Sokoloff, 19 F. Blanc, 20 P. C. Bloom, 20 W. T. Ford, 20 A. Gaz, 20 J. F. Hirschauer, 20 A. Kreisel, 20 

M. Nagel, 20 U. Nauenberg, 20 J. G. Smith, 20 K. A. Ulmer, 20 S. R. Wagner, 20 R. Ayad, 21 'Q A. SofTer, 21 ^ 

W. H. Toki, 21 R. J. Wilson, 21 D. D. Altenburg, 22 E. Feltresi, 22 A. Hauke, 22 H. Jasper, 22 M. Karbach, 22 J. Merkel, 22 

A. Petzold, 22 B. Spaan, 22 K. Wacker, 22 M. J. Kobel, 23 W. F. Mader, 23 R. Nogowski, 23 K. R. Schubert, 23 

R. Schwierz, 23 J. E. Sundermann, 23 A. Volk, 23 D. Bernard, 24 G. R. Bonneaud, 24 E. Latour, 24 Ch. Thiebaux, 24 

M. Verderi, 24 P. J. Clark, 25 W. Gradl, 25 S. Playfer, 25 J. E. Watson, 25 M. Andreotti afc , 26 D. Bettoni a , 26 C. Bozzi a , 26 

R. Calabrese ofc , 26 A. Cecchi ab , 26 G. Cibinetto a6 , 26 P. Franchini a6 , 26 E. Luppi a6 , 26 M. Negrini a6 , 26 A. Petrella afc , 26 

L. Piemontese", 26 V. Santoro a6 , 26 R. Baldini-Ferroli, 27 A. Calcaterra, 27 R. de Sangro, 27 G. Finocchiaro, 27 

S. Pacetti, 27 P. Patteri, 27 I. M. Peruzzi, 27 '! M. Piccolo, 27 M. Rama, 27 A. Zallo, 27 A. Buzzo a , 28 R. Contri a& , 28 

M. Lo Vetere a6 , 28 M. M. Macri a , 28 M. R. Monge ofc , 28 S. Passaggio a , 28 C. Patrignani a6 , 28 E. Robutti a , 28 

A. Santroni a6 , 28 S. Tosi ab , 28 K. S. Chaisanguanthum, 29 M. Morii, 29 R. S. Dubitzky, 30 J. Marks, 30 S. Schenk, 30 

U. Uwer, 30 V. Klose, 31 H. M. Lacker, 31 G. De Nardo a6 , 32 L. Lista Q , 32 D. Monorchio" 6 , 32 G. Onorato afc , 32 

C. Sciacca ab , 32 D. J. Bard, 33 P. D. Dauncey, 33 J. A. Nash, 33 W. Panduro Vazquez, 33 M. Tibbetts, 33 P. K. Behera, ; 

X. Chai, 34 M. J. Charles, 34 U. Mallik, 34 J. Cochran, 35 H. B. Crawley, 35 L. Dong, 35 W. T. Meyer, 35 S. Prell, 35 

E. I. Rosenberg, 35 A. E. Rubin, 35 Y. Y. Gao, 36 A. V. Gritsan, 36 Z. J. Guo, 36 C. K. Lae, 36 A. G. Denig, 37 

M. Fritsch, 37 G. Schott, 37 N. Arnaud, 38 J. Bequilleux, 38 A. D'Orazio, 38 M. Davier, 38 J. Firmino da Costa, 38 

G. Grosdidier, 38 A. Hocker, 38 V. Lepeltier, 38 F. Le Diberder, 38 A. M. Lutz, 38 S. Pruvot, 38 P. Roudeau, 38 

M. H. Schune, 38 J. Serrano, 38 V. Sordini, 38 ^ A. Stocchi, 38 G. Wormser, 38 D. J. Lange, 39 D. M. Wright, 39 

I. Bingham, 40 J. P. Burke, 40 C. A. Chavez, 40 J. R. Fry, 40 E. Gabathuler, 40 R. Garnet, 40 D. E. Hutchcroft, 40 

D. J. Payne, 40 C. Touramanis, 40 A. J. Bevan, 41 K. A. George, 41 F. Di Lodovico, 41 R. Sacco, 41 M. Sigamani, 41 
G. Cowan, 42 H. U. Flaecher, 42 D. A. Hopkins, 42 S. Paramesvaran, 42 F. Salvatore, 42 A. C. Wren, 42 D. N. Brown, 43 

C. L. Davis, 43 K. E. Alwyn, 44 N. R. Barlow, 44 R. J. Barlow, 44 Y. M. Chia, 44 C. L. Edgar, 44 G. D. Lafferty, 44 

T. J. West, 44 J. I. Yi, 44 J. Anderson, 45 C. Chen, 45 A. Jawahery, 45 D. A. Roberts, 45 G. Simi, 45 J. M. Tuggle, 45 

C. Dallapiccola, 46 S. S. Hertzbach, 46 X. Li, 46 E. Salvati, 46 S. Saremi, 46 R. Cowan, 47 D. Dujmic, 47 

P. H. Fisher, 47 K. Koeneke, 47 G. Sciolla, 47 M. Spitznagel, 47 F. Taylor, 47 R. K. Yamamoto, 47 M. Zhao, 47 

S. E. Mclachlin, 48 'Q P. M. Patel, 48 S. H. Robertson, 48 A. Lazzaro afc , 49 V. Lombardo a , 49 F. Palombo a6 , 49 

J. M. Bauer, 50 L. Cremaldi, 50 V. Eschenburg, 50 R. Godang, 5O -0 R. Kroeger, 50 D. A. Sanders, 50 D. J. Summers, 50 

H. W. Zhao, 50 M. Simard, 51 P. Taras, 51 F. B. Viaud, 51 H. Nicholson, 52 M. A. Baak, 53 G. Raven, 53 H. L. Snoek, 53 

C. P. Jessop, 54 K. J. Knoepfel, 54 J. M. LoSecco, 54 W. F. Wang, 54 G. Benelli, 55 L. A. Corwin, 55 K. Honscheid, 55 

H. Kagan, 55 R. Kass, 55 J. P. Morris, 55 A. M. Rahimi, 55 J. J. Regensburger, 55 S. J. Sekula, 55 Q. K. Wong, 55 

N. L. Blount, 56 J. Brau, 56 R. Frey, 56 O. Igonkina, 56 J. A. Kolb, 56 M. Lu, 56 R. Rahmat, 56 N. B. Sinev, 56 D. Strom, 56 

J. Strube, 56 E. Torrence, 56 G. Castelli ah , 57 N. Gagliardi ah , 57 M. Margoni a6 , 57 M. Morandin a , 57 M. Posocco , 57 

M. Rotondo", 57 F. Simonetto afc , 57 R. Stroili ab , 57 C. Voci a6 , 57 P. del Amo Sanchez, 58 E. Ben-Haim, 58 H. Briand, 58 



34 



G. Calderini, 58 J. Chauveau, 58 P. David, 58 L. Del Buono, 58 O. Hamon, 58 Ph. Leruste, 58 J. Ocariz, 58 A. Perez, 58 

J. Prendki, 58 L. Gladney, 59 M. Biasini afc 60 R. Covarelli afc , 60 E. Manoni a6 , 60 C. Angelini afc , 61 G. Batignani afc , 61 

S. Bettarini a6 , 61 M. Carpinelli afc , 61 '[lt] A. Cervelli afc , 61 F. Forti a6 , 61 M. A. Giorgi afc , 61 A. Lusiani ac , 61 

G. Marchiori a6 , 61 M. Morganti ab , 61 N. Neri a6 , 61 E. Paoloni a6 , 61 G. Rizzo a6 , 61 J. J. Walsh Q , 61 J. Biesiada, 62 

D. Lopes Pegna, 62 C. Lu, 62 J. Olsen, 62 A. J. S. Smith, 62 A. V. Telnov, 62 F. Anulli Q , 63 E. Baracchini afc , 63 

G. Cavoto a , 63 D. del Re a6 , 63 E. Di Marco ab , 63 R. Faccini a6 , 63 F. Ferrarotto a , 63 F. Ferroni a6 , 63 M. Gaspero ab , 63 

P. D. Jackson a , 63 L. Li Gioi a , 63 M. A. Mazzoni a , 63 S. Morganti a , 63 G. Piredda a , 63 F. Polci Q& , 63 F. Renga afc , 63 

C. Voena a , 63 M. Ebert, 64 T. Hartmann, 64 H. Schroder, 64 R. Waldi, 64 T. Adye, 65 B. Franek, 65 E. O. Olaiya, 65 

W. Roethel, 65 F. F. Wilson, 65 S. Emery, 66 M. Escalier, 66 L. Esteve, 66 A. Gaidot, 66 S. F. Ganzhur, 66 
G. Hamel de Monchenault, 66 W. Kozanecki, 66 G. Vasseur, 66 Ch. Yeche, 66 M. Zito, 66 X. R. Chen, 67 H. Liu, 67 

W. Park, 67 M. V. Purohit, 67 R. M. White, 67 J. R. Wilson, 67 M. T. Allen, 68 D. Aston, 68 R. Bartoldus, 68 

P. Bechtle, 68 J. F. Benitez, 68 R. Cenci, 68 J. P. Coleman, 68 M. R. Convcry, 68 J. C. Dingfelder, 68 J. Dorfan, 68 

G. P. Dubois-Felsmann, 68 W. Dunwoodie, 68 R. C. Field, 68 A. M. Gabareen, 68 S. J. Gowdy, 68 M. T. Graham, 68 

P. Grenier, 68 C. Hast, 68 W. R. Inncs, 68 J. Kaminski, 68 M. H. Kelsey, 68 H. Kim, 68 P. Kim, 68 M. L. Kocian, 68 

D. W. G. S. Leith, 68 S. Li, 68 B. Lindquist, 68 S. Luitz, 68 V. Luth, 68 H. L. Lynch, 68 D. B. MacFarlane, 68 

H. Marsiske, 68 R. Messner, 68 D. R. Muller, 68 H. Neal, 68 S. Nelson, 68 C. P. O'Grady, 68 I. Ofte, 68 A. Perazzo, 68 

M. Perl, 68 B. N. Ratcliff, 68 A. Roodman, 68 A. A. Salnikov, 68 R. H. Schindler, 68 J. Schwiening, 68 A. Snyder, 68 

D. Su, 68 M. K. Sullivan, 68 K. Suzuki, 68 S. K. Swain, 68 J. M. Thompson, 68 J. Va'vra, 68 A. P. Wagner, 68 

M. Weaver, 68 C. A. West, 68 W. J. Wisniewski, 68 M. Wittgen, 68 D. H. Wright, 68 H. W. Wulsin, 68 A. K. Yarritu, 68 

K. Yi, 68 C. C. Young, 68 V. Ziegler, 68 P. R. Burchat, 69 A. J. Edwards, 69 S. A. Majewski, 69 T. S. Miyashita, 69 

B. A. Petersen, 69 L. Wilden, 69 S. Ahmed, 70 M. S. Alain, 70 R. Bula, 70 J. A. Ernst, 70 B. Pan, 70 M. A. Saeed, 70 

S. B. Zain, 70 S. M. Spanier, 71 B. J. Wogsland, 71 R. Eckmann, 72 J. L. Ritchie, 72 A. M. Ruland, 72 C. J. Schilling, 72 

R. F. Schwitters, 72 B. W. Drummond, 73 J. M. Izen, 73 X. C. Lou, 73 F. Bianchi a6 , 74 D. Gamba ab , 74 M. Pelliccioni afc , 74 

M. Bomben afc , 75 L. Bosisio afc , 75 C. Cartaro Q& , 75 G. Delia Ricca a6 , 75 L. Lanceri afc , 75 L. Vitale Q& , 75 V. Azzolini, 76 

N. Lopez-March, 76 F. Martinez- Vidal, 76 D. A. Milanes, 76 A. Oyanguren, 76 J. Albert, 77 Sw. Banerjee, 77 

B. Bhuyan, 77 H. H. F. Choi, 77 K. Hamano, 77 R. Kowalewski, 77 M. J. Lewczuk, 77 I. M. Nugent, 77 J. M. Roney, 77 

R. J. Sobie, 77 T. J. Gershon, 78 P. F. Harrison, 78 J. Ilic, 78 T. E. Latham, 78 G. B. Mohanty, 78 H. R. Band, 79 

X. Chen, 79 S. Dasu, 79 K. T. Flood, 79 Y. Pan, 79 M. Picrini, 79 R. Prepost, 79 C. O. Vuosalo, 79 and S. L. Wu 79 

(The BABAR Collaboration) 

1 Laboratoire de Physique des Particules, IN2P3/CNRS et Universite de Savoie, F-74941 Annecy-Le-Vieux, France 

2 Universitat de Barcelona, Facultat de Fisica, Departarnent ECM, E-08028 Barcelona, Spain 

8 INFN Sezione di Bari a ; Dipartmento di Fisica, Universita di Barf , 1-70126 Bari, Italy 

* University of Bergen, Institute of Physics, N-5007 Bergen, Norway 

5 Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 

6 University of Birmingham, Birmingham, B15 2TT, United Kingdom 

'Ruhr Universitat Bochum, Institut fur Experimentalphysik 1, D-44780 Bochum, Germany 

8 University of Bristol, Bristol BS8 1TL, United Kingdom 

University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1 

10 Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 

11 Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

12 University of California at Irvine, Irvine, California 92697, USA 

13 University of California at Los Angeles, Los Angeles, California 90024, USA 

1Jh University of California at Riverside, Riverside, California 92521, USA 

15 University of California at San Diego, La Jolla, California 92093, USA 

16 University of California at Santa Barbara, Santa Barbara, California 93106, USA 

17 University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 

18 California Institute of Technology, Pasadena, California 91125, USA 

19 University of Cincinnati, Cincinnati, Ohio 45221, USA 

20 University of Colorado, Boulder, Colorado 80309, USA 

21 Colorado State University, Fort Collins, Colorado 80523, USA 

22 Technische Universitat Dortmund, Fakultat Physik, D-44^21 Dortmund, Germany 

21 Technische Universitat Dresden, Institut fur Kern- und Teilchenphysik, D-0 1062 Dresden, Germany 

2J * Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, F-91128 Palaiseau, France 

25 University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

26 INFN Sezione di Ferrara a ; Dipartimento di Fisica, Universita di Ferrara b , 1-44100 Ferrara, Italy 

27 INFN Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 

28 INFN Sezione di Genova a ; Dipartimento di Fisica, Universita di Genova b , 1-16146 Genova, Italy 

29 Harvard University, Cambridge, Massachusetts 02138, USA 



30 Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

Humboldt- Universitat zu Berlin, Institut fur Physik, Newtonstr. 15, D-12489 Berlin, Germany 

INFN Sezione di Napoli a ; Dipartimento di Scienze Fisiche, 

Universitd di Napoli Federico IF, 1-80126 Napoli, Italy 

Jx Imperial College London, London, SW7 2AZ, United Kingdom 

SJ> University of Iowa, Iowa City, Iowa 52242, USA 

35 Iowa State University, Ames, Iowa 50011-3160, USA 

36 Johns Hopkins University, Baltimore, Maryland 21218, USA 

37 Universitat Karlsruhe, Institut fur Experimentelle Kernphysik, D-76021 Karlsruhe, Germany 

3 Laboratoire de I'Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 

Centre Scientifique d'Orsay, B. P. 34, F-91898 ORSAY Cedex, France 

39 Lawrence Livermore National Laboratory, Livermore, California 94550, USA 

40 University of Liverpool, Liverpool L69 7ZE, United Kingdom 

Jtl Queen Mary, University of London, El 4NS, United Kingdom 

42 University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 

J>3 University of Louisville, Louisville, Kentucky 40292, USA 

^ University of Manchester, Manchester M13 9PL, United Kingdom 

45 University of Maryland, College Park, Maryland 20742, USA 

46 University of Massachusetts, Amherst, Massachusetts 01003, USA 

J>1 Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 

i8 McGill University, Montreal, Quebec, Canada H3A 2T8 

49 INFN Sezione di Milano a ; Dipartimento di Fisica, Universitd di Milano b , 1-20133 Milano, Italy 

50 University of Mississippi, University, Mississippi 38677, USA 

jl Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 3,17 

Mount Holyoke College, South Hadley, Massachusetts 01075, USA 

53 NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

54 University of Notre Dame, Notre Dame, Indiana 46556, USA 

ss Ohio State University, Columbus, Ohio 43210, USA 

56 University of Oregon, Eugene, Oregon 97403, USA 

57 INFN Sezione di Padova a ; Dipartimento di Fisica, Universitd di Padova b , 1-35131 Padova, Italy 

58 Laboratoire de Physique Nucleaire et de Hautes Energies, 

IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, 

Universite Denis Diderot- Paris7, F-75252 Paris, France 

59 University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 

60 INFN Sezione di Perugia a ; Dipartimento di Fisica, Universitd di Perugia b , 1-06100 Perugia, Italy 

61 INFN Sezione di Pisa a ; Dipartimento di Fisica, 

Universitd di Pisa b ; Scuola Normale Superiore di Pisa c , 1-56127 Pisa, Italy 

62 Princeton University, Princeton, New Jersey 08544, USA 

63 INFN Sezione di Roma a ; Dipartimento di Fisica, 

Universitd di Roma La Sapienza b , 1-00185 Roma, Italy 

64 Universitat Rostock, D- 18051 Rostock, Germany 

65 Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

66 DSM/Dapnia, CEA/Saclay, F-91191 Gif-sur-Yvette, France 

67 University of South Carolina, Columbia, South Carolina 29208, USA 

h8 Stanford Linear Accelerator Center, Stanford, California 94309, USA 

69 Stanford University, Stanford, California 94305-4060, USA 

70 State University of New York, Albany, New York 12222, USA 

71 University of Tennessee, Knoxville, Tennessee 37996, USA 

72 University of Texas at Austin, Austin, Texas 78712, USA 
73 University of Texas at Dallas, Richardson, Texas 75083, USA 

7 ^INFN Sezione di Torino a ; Dipartimento di Fisica Sperimentale, Universitd di Torino b , 1-10125 Torino, Italy 

75 INFN Sezione di Trieste" ; Dipartimento di Fisica, Universitd di Trieste b , 1-34127 Trieste, Italy 

76 IFIC, Universitat de Valencia- C SIC, E-46071 Valencia, Spain 

77 University of Victoria, Victoria, British Columbia, Canada V8W 3P6 

78 Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 

79 University of Wisconsin, Madison, Wisconsin 53706, USA 

(Dated: October 3, 2008) 

We report a study of B — > D* K decays with D* decaying to D-k or Dj, using 383 x 10 6 PP> 
pairs collected at the Y(4S) resonance with the BABAR detector at the SLAC PEP-II B-Factory. 
The D meson decays under study include a non-CP mode (K 7r T ), CP-even modes (K K^ , -k 7r T ) 
and CP-odd modes {K®-k°, K®<f>, K%uj). We measure ratios (Rcp±) of branching fractions of decays 
to CP eigenmode states and to flavor-specific states as well as CP asymmetries {Aq P± ). These 
measurements are sensitive to the unitarity triangle angle 7. We obtain Aq P+ — — 0.11±0.09±0.01, 



R* C p+ = 1.31±0.13±0.04, and A* CP _ = 0.06±0.10±0.02, R* CP _ = 1.10±0.12±0.04, where the first 
error is statistical and the second error is systematic. Translating our results into an alternative 
parametrization, widely used for related measurements, we obtain x* + — 0.11 ± 0.06 ± 0.02 and 
x*_ — 0.00 ± 0.06 ± 0.02. No significant CP-violating charge asymmetry is found in either the 
flavor-specific mode D — > K ± n^ 1 or in B — » D*^ decays. 

PACS numbers: 14.40.Nd, 13.25.Hw, 12.15.Hh, 11.30.Er 



I. INTRODUCTION 

In the Standard Model (SM), CP-violating phenom- 
ena are a consequence of a single complex phase in the 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix 0. The B ± -> D^K < -*^ ± decay modes provide 
a theoretically clean determination of the unitarity tri- 
angle angle 7, since the latter is equal to the relative 
phase between the CKM- and color-favored b — * c and 
the CKM- and color-suppressed b — > u decay amplitudes 
that are dominant in the considered decays. The method 
proposed by Gronau, London and Wyler (GLW) makes 
use of the direct CP violation in the interference between 
the amplitudes iov_B ± ->D Q K ± and B ± ^D°K ± decays 
when the D° and D° mesons decay to the same CP eigen- 
state 0, Q- The same approach is equally applicable 
when the D and/or the K meson is replaced with its ex- 
cited state. In this paper we use the B ± — > D*K ± decay. 
We use the notation D°, D*°, D° and D*° to denote 
states with definite flavor, while Dcp+ and Dq P+ de- 
note CP-even eigenstates, Dcp- and B>* CP _ denote CP- 
odd eigenstates, and D and D* denote any state of the 
D(1864)° and £>*(2007)° mesons, respectively. With the 
integrated luminosity presently available, it is not pos- 
sible to make a precise 7 measurement with the GLW 
method alone, but the combination of several methods 
and of several modes allows an improvement of the over- 
all precision ]4|. 

In the case of B^ — ► D*K ± decays, one defines the 
CP-violating charge asymmetry 
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and the ratio of branching fractions for the decays to CP 
eigenmodes and flavor-specific states 
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We refer to the companion of the charmed meson in the 
final state as the prompt track. Experimentally, it is 
convenient to normalize the branching fractions of the 
decays with a prompt kaon in the final state to those 
of the similar decays with a prompt pion that have a 
larger branching fraction 
expressed as 



The ratio R* CP ± can then be 
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where R± and R* are the K/ir ratios 
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predicted to be of the order of 
/V ud \f = 0.080 ±0.002 [|, where f K and 



The ratio R* 

[(fx/U) x \V US 
Jtt are the form factors of the mesons. Equation (J3J) would 
be an equality if CP violation was completely absent in 
B — > D*ir decays. Defining the charge asymmetry 



At = 



B(B- -* D*h~) - B(B+ -> D*h+) 
B{B~ -» D*h~) + B(B+ -> D*h+) 



(6) 



(noted A* and A* K when referring to h = it and h = K 
respectively), this approximation implies that the pion 
charge asymmetry A* should be compatible with zero, as 
should be the kaon charge asymmetry A* K for the flavor- 
specific modes D — > A" ± 7r =F . The possible bias induced 
by this approximation is expected to be small since the 
ratio of the amplitudes of the B~ — > D*°ir~ and B~ — > 
D*°tt~ processes is predicted to be of the order of 1% p] 
in the SM, and will be accounted for in the systematic 
uncertainties. 

Most experimental systematic uncertainties, such as 
those related to the reconstruction of the D* , and the 
uncertainties on the D decay branching fractions, cancel 
in the K/tt ratios R* and R± . By neglecting the small 0, 
8] D°- D° mixing Q and CP violation in D° decays, 
R*cp± and A* CP ± are related to 7 through 
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TABLE I: Past measurements of parameters related to the measurement of 7 in B — > D* K decays by the GLW method. 
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R* 



BABAR [12] -0.10 ±0.23; 

Belle [13, 14] -0.20 ± 0.22 ± 0.04 0.13 ±0.30 ± O.i 



, +0.03 
-0.04 



1.06 ± 0.26 toll 
1.41 ± 0.25 ± 0.06 



0.0813 ± 0.0040 to 0031 
1.15 ± 0.31 ± 0.12 0.078 ± 0.019 ± 0.009 



where r* B is the magnitude of the ratio of the ampli- bias) [211 ]. The parameters x± can be obtained from 



B' 



D* {) K- and B~ 



B is the relative strong phase between 



tudes for the processes 
D*°K-, and S 
these two amplitudes. The ratio r B involves a CKM 
factor IVufjVcis/V^KisI ~ 0.44 ± 0.05 J5| and a color sup- 
pression factor that has been estimated to lie between 
0.26 ± 0.07 ± 0.05 [H and 0.44 Q, so that r* B is pre- 
dicted to be in the range 0.1 — 0.2. More recent calcula- 
tions that take into account final state interactions [ill ] 
yield predictions of r* B = 0.09 ± 0.02. 

The latest results by BaBAR and Belle are reported 
in references [l2| and [1J, LL4( respectively. BABAR used 
123 x 10 6 BB pairs with D* -> Dir° and D reconstructed 
in the CP-even modes K + K~ and ir + ir~, and non-CP 
modes K ± n T , K ± ir^n + TT T and K^n^n . Belle used 
275 x 10 6 BB pairs with D* -> Dtt° and D reconstructed 
in the CP-even modes K + K~ and 7r + 7r~, CP-odd modes 
K°tt°, K° s uj, K° s (j) and non-CP modes K^if? 13]. The 
results are summarized in Table Q] Similar studies have 
been performed on the channels S ± — > DK^ J13l. 1 15L [l6l] 
and B ± -* DK** [l7|. 

The BABAR [l8| and Belle [lj| experiments have re- 
cently obtained estimates of r B and 8 B parameters from 
the overlap of the D° and D° decays in the Dalitz planes 
of some three-body D decays. BABAR obtains r* 



0.135 ± 0.051 and 6% 



-63 



+28n 
30 ) 



'B — 

while Belle obtains 



0.21±0.08±0.02±0.05and^ = (342 



+21 



23 ±4±23)° 

(where the first error is statistical, the second is the ex- 
perimental systematic uncertainty and the third reflects 
the uncertainty on the D decay Dalitz models). 

In this paper, by using (383 ± 4) x 10 6 BB pairs, we 
update the results of our previous study of B ± ^D*K ± 
decays [12[ for D decays to the CP-even modes K + K~, 
7r + 7r~ and to the flavor-specific modes K ± ir^ z , and we ex- 
tend it to the CP-odd modes K®ir°, K®lu and K®(f), and 
to D* — ► D^f. Due to parity and angular-momentum con- 
servation, the CP eigenvalue of the D* is inferred from 
that of the D. when the CP eigenvalue of the neutral com- 
panion (7 or 7T°) is taken into account [20]: CP(D*) — 
CP(D) when D* -v Dn°, and CP(D*) = -CP(D) when 
D* -> Z>y. 

Experimental results can also be presented using the 
"cartesian coordinates" 

(x*±,y±) = (r B cos(S* B ± j),r* B sm(5* B ± 7)), (9) 

which have the advantage of having Gaussian uncertain- 
ties, and of being uncorrelated and unbiased (r B , be- 
ing positive, is biased towards larger values in low preci- 
sion measurements, whereas x J. and y± show no such 



R-cp± ano - A CP ±, 



R 



CP^ 



(IT A 



CP^ 



R* CP _(1 T A* CP _) 



(10) 



The measurements presented in this paper have no direct 
sensitivity to y±, in contrast to Dalitz analyses. However 
an indirect constraint can be obtained using 



r * \2 *2 1 *2 

(xb) = x ± +y± = 



p. 



CP^ 



P, 



CP- 



(11) 



Note that there are four observables in these parameter- 
izations, either (A* CP+ , R* CP+ , A* CP _ and R* CP _) or {x* + , 
y*L, x*_ and y*_), while there are only three independent 
fundamental parameters (7, r B and S B ). The set of pa- 
rameters must therefore fulfill one constraint, which can 
be k = 0, where 



k — Rcp+^cpa 



A* /?* 
j± cp _jx cp _. 



(12) 



II. THE DATASET AND BABAR DETECTOR 

The results presented in this paper are based on 
data collected with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ storage ring of the Stanford 
Linear Accelerator Center. At PEP-II, 9.0 GeV elec- 
trons and 3.1 GeV positrons collide at a center-of-mass 
energy of 10.58 GeV, which corresponds to the mass of 
the T(AS) resonance. The asymmetric beam energies re- 
sult in a boost from the laboratory to the center-of-mass 
frame of (3j ~ 0.56. The dataset analyzed in this paper 
corresponds to an integrated luminosity of 347 fb _ at 
the T(4S) resonance. 

The BABAR detector is described in detail else- 
where [22]. Surrounding the interaction point is a five- 
layer double-sided silicon vertex tracker (SVT), which 
measures the trajectories of charged particles. A 40- 
layer drift chamber (DCH) provides measurements of 
the momenta of charged particles. Both the SVT and 
DCH are located inside a 1.5 T magnetic field provided 
by a solenoid magnet. Charged hadron identification is 
achieved through measurements of particle energy-loss in 
the tracking system and the Cherenkov angle obtained 
from a detector of internally reflected Cherenkov light 
(DIRC). A CsI(Tl) electromagnetic calorimeter (EMC) 
provides photon detection, electron identification, and 
7T° reconstruction. Finally, the instrumented flux return 



(IFR) of the magnet enables discrimination of muons 
from pions. For the most recent 134 fb _1 of data, a frac- 
tion of the resistive plate chambers constituting the muon 
system has been replaced by limited streamer tubes [23| . 
We use Monte Carlo (MC) simulation to study the 
detector acceptance and backgrounds. The MC uses the 
EvtGen generator [H and GEANT4 (2^ to simulate the 
passage of particles through matter. 



III. RECONSTRUCTION OF B CANDIDATES 

We perform an exclusive reconstruction of the full B 
meson decay chain, in the modes described in the intro- 
duction, starting from the final stable products (charged- 
particle tracks and neutral electromagnetic deposits in 
the EMC). 

The 7T° candidates used to form an w, a I? or a D* can- 
didate are reconstructed from pairs of photons with en- 
ergies larger than 30 MeV, and shower shapes consistent 
with electromagnetic showers, with invariant mass in the 
range 115 < m 77 < 150MeV/c 2 . In addition, the tt° can- 
didates used to form a D* candidate are required to have 
center-of-mass frame momenta p* < 450MeV/c. The 
u) candidates are reconstructed in the three-body decay 
u> — > 7r + 7r~7r , with an invariant mass within 50MeV/c 2 
of the world average Q. We reconstruct K® — > 7r + 7r~ 
candidates from pairs of oppositely charged tracks that 
are consistent with having originated from a common ver- 
tex position and with an invariant mass within 25 MeV/c 2 
of the world average Q. We reconstruct (f> — ► K + K~ 
candidates from pairs of oppositely charged tracks with 
particle identification inconsistent with a pion hypoth- 
esis, that are consistent with having originated from a 
common vertex position, and that have invariant mass 
within 30 MeV/c 2 of the world average Q • 

Only two-body D decays are considered in this study. 
The D candidates are reconstructed from their two 
daughters that are required to be consistent with hav- 
ing originated from a common vertex position. In the 
case of D — ► K®tt , in which vertexing of the K®tt° sys- 
tem would yield a poor geometrical constraint, a beam 
spot constraint is added to the fit in order to force the D 
daughters to originate from the interaction region. 

The D* candidates are formed from D and tt° or 7 
candidates. These photon candidates are required to 
have energies larger than 100 MeV and shower shapes 
consistent with electromagnetic showers. The D* can- 
didates are required to fulfill 130 < Am < 170MeV/c 2 
and 80 < Am < 180MeV/c 2 , respectively, where Am is 
the invariant mass difference between the D* and the D 
candidate. 

The 7T°, K®, D and D* candidates are refitted with 
mass constraints before their four-momenta are used 
to reconstruct the B decay chain. We form B can- 
didates from D* candidates and charged tracks, fit- 
ted with a beam spot constraint. We characterize 
B candidates by two kinematic variables: the differ- 



ence between the reconstructed energy of the B candi- 
date and the beam energy in the center-of-mass frame 
A Ek = Eg — y/s/2, and the beam-energy substituted 
mass m-Es = y/(s/2 + p • p B ) 2 /£$ - p|, where (E , 
Po) and (Eb, Pb) are the four-momenta of the T(AS) 
and B meson candidate, respectively, the asterisk de- 
notes the T(4S f ) rest frame, and -^/s is the total energy 
in the T(4S) rest frame. The subscript K in AEk in- 
dicates that the kaon hypothesis has been assumed for 
the prompt track in the computation of AE. For a cor- 
rectly reconstructed B meson having decayed to a D* K 
final state, AEk is expected to peak near zero, with an 
R.M.S. of about 16 MeV, and toes is expected to peak 
near the B-meson mass 5.279 GeV/c 2 , with an R.M.S. 
that is almost independent of the channel and close to 
3 MeV/c 2 . For a B — » D*ir decay reconstructed as 
B — > D*K with a correctly identified D* , the AE K peak 
is shifted by approximately +50 MeV. At reconstruction 
level, the loose requirements 5.2 < toes < 5.3 GeV/c 2 and 
\AEk\ < 0.2 GeV are applied to the B meson candidate. 
We form a Fisher discriminant T [26] to distinguish 
signal events from the significant background due to 
e + e~ — ► qq (q = u,d,s,c) continuum events. Six vari- 
ables are used: 

• Lq and Li , the zeroth and second angular moments 
of the energy flow around the B thrust axis. They 
are defined as YliPi an d JVpiCOS 2 0i respectively, 
where pi is the momentum and 9i is the angle with 
respect to the thrust axis of the B candidate, both 
in the center-of-mass frame, for all tracks and neu- 
tral clusters not used to reconstruct the B meson; 

• R2, the ratio of the second to the zeroth Fox- 
Wolfram moment [27| of charged tracks and neutral 
clusters in the center-of-mass frame; 

• I cos 9b I , where 0b is the angle between the momen- 
tum of the B candidate and the boost direction of 
the e + e~ center-of-mass frame; 

• J cos ^Thrust 1 1 where ^Thrust is the angle between the 
B candidate thrust vector and the beam axis in the 
center-of-mass frame; 

• I cos #t I j where #t is the angle between the B can- 
didate thrust axis and the thrust axis of the rest 
of the event in the center-of-mass frame (where the 
rest of the event corresponds to reconstructed par- 
ticles not associated with the B candidate). 



IV. SELECTION OF B CANDIDATES 

After the preliminary event reconstruction, a large 
amount of background remains in the signal candidate 
sample. In this section we describe the additional selec- 
tion criteria used to reduce the background. 

The selection of each B — > D* K decay mode is opti- 
mized separately, by the maximization of the sensitivity 



S/y/S + B + 1, where V5 1 + 73 + 1 is a symmetrized ap- 
proximation of the Poisson uncertainty on the measure- 
ment of S + B. The numbers S and B of signal and back- 
ground expected events are estimated from, respectively, 
high-statistics exclusive MC samples, and a cocktail of 
generic B + B~ (with signal events removed), B°B° and 
qq MC samples. 

In the optimization procedure, we include require- 
ments on all variables, including those that will be re- 
laxed during the fit, and including tightening require- 
ments that have been made in the reconstruction stage. 
Our optimization procedure is similar to that used in 
Ref. [28|, and allows us to determine the optimal set of 
variables as well as the optimal requirements on those 
variables, by the examination of the signal-to-background 
ratio distributions [29| , The final set of variables on 
which we apply selection optimization is: 

• the B candidate-related variables jjjes an d AEk 
introduced above; 

• the mass m D o of the D candidate before the mass 
constraint is applied and the mass difference Am; 

• likelihood ratios for the prompt track, that are eval- 
uated making use of the Cherenkov angle infor- 
mation from the DIRC, and of the dE/dx infor- 
mation provided by the tracking system. Explic- 
itly, we compute likelihoods £h for particle identi- 
fication (PID) hypotheses h — K,ir,p and make 
requirements on the ratios £k/{£k + £■*) and 
£k/{£k + £ P )- We also require that the track is 
not identified as an electron or muon; 

• likelihood ratios for pion and kaon candidates that 
are daughters of two-body D decays; 

• the value of the Fisher variable T\ 

• the invariant masses of the K®, </>, 7r° and u> can- 
didates, when relevant, and before the mass con- 
straints. Furthermore, for decays involving K® can- 
didates, we include the ratio of the flight length 
of Kg candidates in the transverse plane divided 
by its uncertainty, and require it to be larger than 
two. For decays involving u> candidates, we include 
| cos(# w )|, where 9^ is the angle between the normal 
to the pion decay plane and the D direction in the 
to rest-frame. 

The selection requirements applied to these variables are 
mode-dependent, except for the prompt track PID re- 
quirements £k/{£k + £-n) > 0.9 and £k/(£k + £p) > 
0.2 that are applied to all B — ► D*K channels. 
The selection of the B^ — > D*n ± modes is identical to 
that of the i? ± — > D*K ± modes, except for the prompt- 
track PID that is reversed {£ K /(£ K + £„) < 0.2). 

A fraction of the events have several B candidates se- 
lected: the average multiplicity varies from 1.07 to 1.66 
for D* -> Dtt° and from 1.00 to 1.25 for D* -> 7J 7 , de- 
pending on the channel. We select the B candidate that 



has the B vertex fit with the largest probability. This 
best-candidate procedure is used during the optimiza- 
tion of the selection that we have described above. The 
probability of selection of the well-reconstructed signal 
candidate is mode dependent and is in the range 56- 
72% for D* -> L»vr° decays and in the range 68-81% for 
D* — > D-f decays. 



V. MAXIMUM LIKELIHOOD FIT 

The dominant contribution to the remaining back- 
ground after event selection is from B decays, including 
a significant amount of feed-across from B^ — > D*n ± de- 
cays. Therefore the measurement is performed with an 
unbinned likelihood fit [30, [3l[ based on two variables 
that best discriminate this background, namely AEk 
and a PID variable T-ji defined below. 

As the prompt track PID likelihood ratio 72 = 
£k/(£k + £tt) is very strongly peaked near zero for pi- 
ons and near one for kaons, we use a pseudo-logarithmic 
change of variable 



Tn = log 



10 



K 



l-K 



(13) 



We include a small positive number e = 10 , so that 
T-ji is always defined, with 7^ = +7 for TZ = 1 ("perfect 
kaons") and %^ = — 7 for 71 = ("perfect pions"). 

The measurement is performed with an extended un- 
binned maximum likelihood function 



-JV' 



£ 



(AT 



N\ 






(14) 



8=1 



where N is the number of events in the sample to fit, N' 
is the expected number, and for event i 



* = ^£^. 



(15) 



where j = D*K, D*ir. Bk, B n is one of four events cate- 
gories: signal kaon and pion, and background kaon and 
pion respectively, where the background is a combination 
of continuum, B + B~ and B°B events. The quantity V\ 
is the probability density function (PDF) for event i and 
category j, and Nj is the number of events in category j. 

For the signal categories, the distance between the 
kaon and pion AEk peaks provides powerful separa- 
tion between pions and kaons, in addition to PID. For 
the background categories, we use mutually exclusive 
likelihood-based pion and kaon selectors, that in particu- 
lar contain requirements of 72. > 0.9 (kaon) and 72. < 0.1 
(pion) respectively. For consistency and symmetry rea- 
sons, the whole region 0.1 < 72 < 0.9 is removed for all 
categories, including the signal categories used in the fit. 

The correlations between T-ji and AE^ are found to 
be small (compatible with zero for the signal K and for 




-0.05 0.05 



0.1 -0.05 0.05 0.1 




AEk(GoV) 




210 J r 



-6-4-2024 




■0.05 0.05 0.1 



i 

I io-' 



1 10-' 



-20246 



Tn 

FIG. 1: Distributions of AEk (upper plots) and the PID variable Tn (lower plots) from MC simulations of the categories (from 
left to right) D*K, D*<k, B K and B n (the latter two from B + B~ , B°B° and qq MC samples), for the mode D* -» D%°,D -» 
K 7r T . In the upper plots, the dots represent the MC sample spectrum, and the curves show the PDFs. Note the vertical log 
scale in the lower plots. 



the background categories, and with —6% for the 7r signal 
category) therefore a factorized approximation is used 



V\ (AE K , T n ) = V\ (AE K )Vf (T n ). 



(16) 



We have checked that no bias is introduced by this ap- 
proximation by simulating a large number of experiments 
in which the signal is taken from the large statistics exclu- 
sive MC samples used for estimating these correlations. 

The PDFs used in the fit are determined from MC 
samples. The signal AEk PDFs are parameterized with 
double Gaussian functions. The background PDFs are 
mode-dependent functions chosen to best represent the 
MC background distributions: they include Gaussian, ex- 
ponential, and third-order Chebyshev polynomial func- 
tions. The complicated shape of the AEk distribu- 
tion of the B^ category arises from the contributions 
from several distinct components: at low AEk values, 
S ± — ► D* p^ decays dominate; in the signal region, the 
background is mainly composed of 7 <-> 7r° cross- feed and 
of B ± — » D^n^, the latter of which dominates at high 
AEk values. The T-R PDFs are histograms, determined 
from MC samples, with a binning AT-ji — 0.5 and, there- 
fore, 28 bins. MC-based studies have shown that the 
results of such fits do not depend on the number of bins 
rib as long as n& > 2. 

We correct for a small discrepancy in PID efficiencies 
between data and MC, using high-statistics high-purity 
kaon and pion samples from inclusive D*^ — > Dtt , D — v 
K ± tt^ data. The difference in track momentum spectra 
between these control samples and the exclusive modes 
studied in the present analysis is accounted for in the 
correction procedure. This is achieved by weighting the 
control sample T-ji PDF by the ratio of the MC to control 
sample prompt track momentum distributions for both 
cases of the prompt track being a kaon or a pion. An 



example of the PDFs used for the channel D* — > Dir , 
D — ► K ± tt t is shown in Fig. [1] 

For signal events with a pion prompt track, for which 
AE n (the subscript n indicates that the pion hypothesis 
has been assumed for the prompt track in the computa- 
tion of AE) is close to zero, the relation 



AE 



K 



AE, 



1 E 7 



2p m r(4S) 



(ml-ml), (17) 



introduces a mild dependence of AEk on the momen- 



tum p of the prompt track. The parameters E- 



T(4S) 



and 



TO r(4S)j TriKi ttI'tt denote the energy of the e + e~ system 
in the laboratory frame and the masses of the mesons, re- 
spectively. Fits taking this dependence into account do 
not show any significant improvement, nor degradation. 

Fits performed on the B + B~ , B°B° and qq background 
MC samples show no significant bias. Similar fits with ei- 
ther pion or kaon signal events removed yield numbers of 
signal events compatible with zero for the removed cate- 
gory. This indicates that the factorization approximation 
made for the background PDF does not create any bias 
on the number of fitted signal events. 

Signal efficiencies are estimated from fits on high statis- 
tics exclusive MC samples and summarized in Table ITT] 
We perform separate fits for each D* decay mode, and 
subsequently perform a weighted average to obtain our 
final results for Rc P ± and Aq P: j_ . The free parameters of 
each fit are itemized here: 

• the charge-averaged Kjix ratio (one parameter, R± 
or R*, whenever relevant); 

• the number of pion signal events (one parameter); 

• the pion and kaon charge asymmetries A* h of the 
signal (two parameters); 



TABLE II: Selection efficiencies (in %) for channels used in 
this analysis for each decay mode of the D (statistical uncer- 
tainties only). 





(Dk )^ 


(D^K* 


(Dk )^ 


(D-y)** (in %) 


#±71"* 


21.0 ± 0.1 


24.7 ± 0.1 


22.2 ± 0.1 


24.9 ± 0.1 


KK 


14.6 ± 0.1 


14.7 ± 0.1 


14.8 ± 0.1 


14.8 ± 0.1 


KK 


20.4 ± 0.1 


21.1 ± 0.1 


20.5 ± 0.1 


21.2 ± 0.1 


KW 


8.9 ± 0.1 


8.8 ± 0.1 


8.9 ± 0.1 


9.0 ± 0.1 


K s lo 


4.4 ± 0.1 


4.2 ± 0.1 


4.5 ± 0.1 


4.3 ± 0.1 


KU 


10.3 ± 0.1 


13.5 ± 0.1 


10.4 ± 0.1 


13.7 ± 0.1 



• the number of pion and kaon background events, 
and charge asymmetries (four parameters); 

• the position of the AEk peak of the pion events 
(one parameter). 

In total there are nine free parameters for each D* mode. 



VI. SYSTEMATIC UNCERTAINTIES 



The systematic uncertainties are summarized in Table 



The contribution of the determination of the AEk sig- 
nal PDFs to the systematic uncertainty is estimated by 
varying the parameters that are fixed in the fit by one 
standard deviation (±ltr). The contribution of the AEk 
PDFs of the B K category for the D*° — » Dir° decays is 
determined similarly. 

For the other AEk background PDFs, this approach 
cannot be used, as the correlations between parameters 
are not small. The contribution to the systematics due to 
the limited MC statistics used to determine the param- 
eters of the PDFs is obtained in the following way. We 
determine two separate parameterizations of the PDFs 
on two halves of the MC sample, and perform the fits 
with them. We take half the difference between the ob- 
tained results as an estimate of the systematics. 

The contribution of the determination of the T-n sig- 
nal PDFs to the systematic uncertainty is estimated by 
performing an additional fit without the correction of the 
small discrepancy between data and MC described above. 
The difference between the results of both fits is taken as 
an estimate of the uncertainty. 

The systematic uncertainty introduced by the limited 
knowledge of B ± — > D* p ± and B° — > D* + ir~ branching 
fractions is estimated from MC samples by performing a 
fit on a sample in which the number of these events is 
varied by ±1<7 [a]. 

Differences in the interactions of positively and nega- 
tively charged kaons with the detector and the possible 
charge asymmetry of the detector are studied using the 
exclusive MC samples. Asymmetries of (—1.0 ± 0.2)% 
and (0.2 ± 0.2)% are observed for kaon and pion modes, 
respectively, for the CP modes. A correction of +1% is 



TABLE III: Contributions to systematic uncertainties for 
each mode on the measurement of the charge asymmetries 
A* K , and the ratio R"cp of CP eigenmode to flavor specific 
mode (10 -3 ). See text for details. 



D* 



Dk u 



K° s uj K° s 4> (10~ 



AE K 





22 


5 


9 


6 


10 


Tn 


1 


12 


2 


18 


22 


38 


B(D*p-) 


3 


4 





1 


2 


59 


B(D*-n+) 


1 


1 


1 


1 


6 





n ^7 





5 


5 


5 


5 


5 


S-wave 


- 


- 


- 


- 


33 


2 



Total 



3 26 7 



21 41 



71 



R*c 



Kk kk KK K° s k° K° s lo K° s 4> (10 -3 ) 



AE K 


80 


34 


54 


116 


54 


T n 


16 


13 


10 


32 


6 


B(D*p-) 


51 


22 


14 


5 


131 


B{D*~k+) 


3 


5 


1 


21 


4 


S-wave 


- 


- 


- 


141 


40 



Total 


96 42 57 187 147 


D* -> D-y 


A K 


K-k kk KK K° s k° K° s uj K° s (j> (10~ 3 ) 



AE K 


1 


39 


8 


19 


66 


52 


Tn 


16 


4 


39 


18 


75 


22 


B(D*p~) 


4 


1 





1 


18 


1 


B{D*-tv+) 


4 


3 


4 





18 


3 


7T <-» 7 





10 


10 


10 


10 


10 


S-wave 


- 


- 


- 


- 


40 


5 



Total 



17 40 41 



28 111 



■ r »s 



B, 



K-k kk KK K%k° K° s oj K° s (j> (10~ 3 ) 



AE K 


- 136 


59 


64 


393 


230 


Tn 


17 


34 


4 


78 


42 


B(D'p-) 


1 


2 


12 


5 


23 


B(D*-k+) 


11 


9 


3 


13 


8 


S-wave 


- 


- 


- 


192 


47 



Total 



138 69 



65 445 



239 



applied to the measured values of A* CP . The simulation 
of the detector charge asymmetry has been compared to 
the actual value in the data in a previous analysis of B 
decays to K-k 32]. The possible discrepancy has been 
found to be smaller than 1%. 

The 7T° <-> 7 cross-feed can reduce the value of A* CP j_ 
because for a given Dqp final state, D^ p has the same 
CP value as Dqp if decaying to Dn° and the opposite 
CP value if decaying to D-y ;2Q]. This "CP dilution" is 
estimated from MC samples by performing a fit in which 
the potential feed-across has been completely removed. 
The effect is similar among modes and is accounted for 
by an uncertainty of 0.5% for Dir° modes and of 1.0% 
for D-y modes 



In the case of D decays to K^ 



and K^u, the CP- 



10 



TABLE IV: Summary of measurements of the charge asymmetries A* K ; the CP eigenmode to flavor specific mode ratios Rqp, 
and the cartesian parameters x* , for B decays to CP even and CP odd eigenmodes D CP K . The value of A* K for the 
flavor-specific control mode with D — » K if* is also given. 





A* K 


Rcp x 


Flavor specific 

CP+ 

CP- 


-0.06 ±0.04 ±0.01 

-0.11 ±0.09 ±0.01 

0.06 ±0.10 ±0.02 


1.31 ±0.13 ±0.04 0.11 ±0.06 ±0.02 
1.10 ±0.12 ±0.04 0.00 ±0.06 ±0.02 



violating charge asymmetry can be diluted by the pres- 
ence of decays to the same final state that may have a 
different CP composition (K®K + K~ and K®ir + ir~ir , re- 
spectively). This S-wave effect is accounted for in a way 
similar to that used in our previous study of the DK 
modes [la ]. It consists of applying a correction to the 
measured A* CP ± and Rqp± values using the CP content 
information of K®K + K~ and K®tt + tt~ tt° modes. The 
uncertainty on the correcting factors is then propagated 
to the correction formula and included as an additional 
systematic uncertainty on A* CP ± and Rg P ±. 

The correlations between the different sources of sys- 
tematic errors are negligible and neglected when combin- 
ing the two CP-even or the three CP-odd modes. 

No systematic error or correction is applied to account 
for selection efficiency uncertainties as we do not measure 
branching fractions but ratios of branching fractions in 
which they largely cancel. 

For the branching fraction ratios R* CP ±, in addition 
to the sources of systematic uncertainties listed in Ta- 
ble IIIII we associate one more uncertainty with the as- 
sumption that Rq P ± — R±/R*. This assumption holds 
only if the magnitude of the ratio r* between the am- 
plitudes of the B~ -> D*°tt- and B~ -> D*°n- pro- 
cesses is neglected j6|. The ratio r* is expected to be 
small: r* ~ r B ^ x2 , where A w 0.22 [a] is the sine of the 
Cabibbo angle. This introduces a relative uncertainty of 
±2r* cos 6* cos 7 on Rq P ± , where 6* is the relative strong 
phase between A(B~ — ► D* Q ir~) and A(B~ — > L>*°7r~). 
Since | cos J* cos 7I < 1 and r* < 0.007, we assign a rela- 
tive uncertainty of ±1.4% to Rc P ±, which is completely 
anti-correlated between R* CP+ and R*q P _ ■ 



A* cp _ n = 0.032 ± 0.027 ± 0.006, and for flavor-specific 
modes At = -0.004 ± 0.010 ± 0.001. 



TABLE V: Number of events measured in this analysis (sta- 
tistical uncertainties only). 





(D-k )^ 


(D^K* 


(Dn )^ 


(£)7) 7 r ± 


JfVF 


874 ± 44 


536 ± 36 


10729 ± 133 


7238 ±119 


7T7T 


31 ±8 


15 ±6 


262 ± 20 


170 ± 17 


KK 


101 ±14 


62 ±12 


987 ± 43 


709 ± 37 


K%^ 


86 ±14 


62 ±11 


900 ± 38 


583 ± 33 


K u s u> 


43 ±11 


29 ±9 


419 ±31 


250 ± 24 


kU 


19 ±6 


21 ±6 


262 ± 22 


180 ± 20 



We also obtain R* = 0.0802 ± 0.0031 ± 0.0018, com- 
patible with the theoretical prediction given in the in- 
troduction, and {r* B ) 2 = 0.20 ± 0.09 ± 0.03. We find 
K = -0.08 ± 0.16 ± 0.02 (defined in Eq. H2J), consistent 
with zero as expected. We confirm the large value of 
(r B ) 2 that can be inferred from the previous measure- 
ments [13| based on the GLW method, with a precision 
improved by a factor two. 

Using the value of 7 = 67.6 ± 4.0 obtained by a SM- 
based fit of the CKM matrix [33| and the values of r B 
and S B from Ref. [ij|, we predict A* CP+ = -0.18 ± 



0.10, R 



CP+ 



1.06 ± 0.06, A 



CP- 



0.20 ± 0.10 and 



R* CP _ = 0.98 ± 0.05. Our results are compatible with 
these predictions. 

We also compute the cartesian coordinates with the 

channel D — > K®4> removed, so as to facilitate the com- 

arison with results obtained with the Dalitz method 



VII. RESULTS 

We plot the AEk distributions in Fig. [2] with a kaon 
selection (7^. > 0) applied, and with the fitted PDFs 
overlaid. The results are summarized in Table IIV( with 
the observed numbers of charged-averaged events in Ta- 
ble [V] Note that none of the corrections between data 
and MC that would be needed for measurements of ab- 
solute branching fractions are used here, as we are in- 
terested in ratios only. We have checked that charge 
asymmetries of B — » D*ir modes are compatible with 



(18) 



zero as expected: A 



CP+,7T 



= 0.007 ± 0.029 ± 0.005, 



x* + = 0.09 ±0.07 ±0.02, 

x*_ = -0.02 ±0.06 ±0.02, 

{r B f = 0.22 ±0.09 ±0.03. 

VIII. SUMMARY 



We have performed measurements of the CP eigen- 
mode to flavor specific mode ratios and of the CP- 
violating charge asymmetries of S ± — > D*K ± decays. 
The ratios Rq P are found to be compatible with, and 
more precise than, previous measurements. Our results 
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FIG. 2: AEk distributions, with a cut on the PID variable Tr > to enhance the kaon part of the sample. Distributions are 
shown for each of the decays modes: (left) D* — > Dtt° and (right) D* — > Dy, with the D decay modes indicated on the left of 
the figure. Dots denote the distribution of the data. Curves denote the PDFs for the various categories: signal K (long-dashed 
curve); signal -n (dotted curve); background K (short-dashed curve); background 7r does not appear because the Tp, > cut 
completely removes this category. The thick curve denotes the total PDF. 
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for Rqpj- and A* CP± are at least a factor of two more 
precise than previous measurements [12, ll3| ■ The preci- 
sion of our results for x± is comparable to that obtained 
from Dalitz plot analyses [Tj| [l9| • No significant charge 
asymmetry is observed in the pion modes. These results 
supersede our previous measurements [12l j . 
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